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Abstract 

We introduce an approach, based on the coordinate space Faddeev 
equations, to solve the quantum mechanical three-body Coulomb problem 
in the continuum. We apply the approach to compute measured properties 
of the first two 0"*" levels in ^^C, including the width of the excited level, 
in a 3a-particle model. We use the two-body interaction which reproduces 
the low energy a-a scattering data and add a three-body force to account 
for the compound state admixture. 



Introduction. The quantum mechanical three-body problem with Coulomb 
interaction is yet one of the unsolved problems of few body physics (see, for 
example, fl], H, |^, ^). A number of three-body Coulomb problems in nuclear, 
atomic and molecular physics cannot be dealt with at the moment because of 
the principal difficulties in constructing correct asymptotic wave functions in the 
case when none of the two body subsystems is bound. 

One of such problems is the description of the astrophysically significant sec- 
ond 0+ state in ^^C (denoted 02") as a system of three alpha particles 0, |^. 
Despite the long history of investigations it is yet unclear to what extent the 
three-body picture accounts for the real structure of this state as it lies above 
the threshold and requires therefore correct continuum spectrum solutions of the 
three-body Coulomb problem. 
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The purpose of this letter is two-fold. First we introduce an approach to 
compute the three-body continuum spectrum for Coulomb interactions. We con- 
struct the asymptotic functions employing the adiabatic hyperspherical expansion 
of the coordinate space Faddeev equations 0]. This method proved to be very 
powerful in treating the long-distance behavior of three-body systems including 
the pathological Efimov effect [Q. The approach is a generalization of the adi- 
abatic hyperspherical method, successfully used in atomic physics @, into the 
more versatile Faddeev equations. 

Secondly, we apply the invented method to the two lowest 0+ states in ^^C 
where descriptions in terms of three interacting a-particles traditionally are be- 
lieved to be rather accurate. The second state is unbound with respect to 
three free a-particles and its structure (along with the ground state structure) 



has been debated for many years, see for example |T^, [IT], 0. A correct descrip- 
tion, which must include continuum effects of the long-range Coulomb interaction, 
is therefore now possible. 

As repeatedly suggested we shall include a genuine three-body force in order 
to account for three-body effects related to the Pauli principle and polarization 
effects [0. Beside providing information about the structure, the method also 
enables us to give an estimate of the width of the 0^-level within the genuine 
three-body model. 

Method. We solve the Faddeev equations in coordinate space using the adi- 
abatic hyperspherical expansion method 0, We shall use hyperspherical 



coordinates, which consist of one radial coordinate p (hyperradius) and five 
generalized angles denoted by O. The hyperradius p is defined in terms of 
the center of mass coordinates of the a-particles as p^ = AaJ2i=irf, where 
Aa = nia/m = 3.97 is given by the nucleon mass m and the a-particle mass ma- 
For a fixed p we consider the eigenvalue equation for the five dimensional 
spin-angular part of the Faddeev operator, i.e. 

(Tn-A„)$51) + |?pV.$r^ = ,^ = 1,2,3 , = , (1) 

where Tq, is the angular part of the kinetic energy operator and Vi is the potential 
between particles j and k, where {i,j,k} is a cyclic permutation of {1,2,3}. 
The eigenvalues A„(p) and the eigenfunctions $;^*^(p, fi), which also carry the 
dependence of the spins of the particles, are calculated as function of p. The total 
three-body wave function is then expanded in terms of these eigenfunctions, i.e. 

P n 

The system of coupled hyperradial equations for the expansion coefficients 
fn{p) finally becomes 

( - 1^ - ^ + 7^ + ^)fniP) = En' {^Pnn'i + Qnn')fn'{p) , (3) 
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For large p the non-diagonal terms P and Q are small compared to the diagonal 
terms A„(p) which, divided by p^, serve as an effective adiabatic hyperradial 
potential. 

The large distance tail of the eigenvalues A„(p) is independent of the properties 
of the short-range nuclear force. The derivation of this asymptotic behavior is 



difficult and beyond the scope of this letter, but the result can be written as [15 



A„(p ^ oo) ^ 3?7p + C„^/?/p + ... , (4) 



where rj = ^{2eY\^A^, 2e is the charge of the a-particle and Cn are constants. 
For the lowest term this constant is Ci = 9/4. At large distances the adiabatic 
potential A„/ p^ has therefore a Coulomb term 3t]/ p, with the next term behaving 
as p~^/^. These terms are absent for uncharged particles. 

For a positive energy E = Ti^K^/{2m) the general solution of equation (j^) 
asymptotically is a linear combination of two linearly independent solutions exp(±i(/tp — |^ log Kp)). 
A resonance with energy Er and width F corresponds to a solution of the system 

with the complex energy Eq = E^ — iV /2 This solution satisfies the 

boundary condition 

fn{p -> cx)) -> exp(+i(/top - — log Kop)), (5) 



where = y2m{Er — iV /2) /fi' . This boundary condition determines that the 
S"- matrix has a pole at the complex energy Eq. 



Interactions. We use the Ali-Bodmer potential version "a" slightly altered 
to reproduce the s-wave resonance of ^Be at 0.093 MeV in addition to the low 
energy aa phase shifts. The potential is given as 

V = (l25Pi=o + 20P/=2)e-'^'/'-''' - 30.18 e""'/'-^^', (6) 



where lengths and energies are in units of fm and MeV, respectively. The P/ is 
the projection operator onto the state with relative orbital angular momentum /. 



It is well known that this potential underbinds the O"*" states |10] due to the 
strong repulsive core which is necessary to simulate the Pauli repulsion in the 
2a system. The width of the resonance depends crucially on its energy, and 
to estimate this width accurately we need the resonance energy at the correct 
position. We achieve this by indtroducng an additional attraction. 

Indeed, when three a-particles are close together one of them polarizes the 
other two beyond the polarization included in the two-body interaction. Thus 
the Pauli repulsion between them should be reduced. This effect can be imitated 
by a phenomenological three-body attractive force, which must be non-vanishing 
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only when all three particles are close together, that is for small p. This "off 
shell" effect must disappear at large p where the "on shell" properties of the 
phase fitted interactions provide the correct asymptotic behavior. As a novel 
feature we include this additional interaction by parametrizing the three-body 
potential as a simple gaussian, V3(p) = 6*3 exp(— which is added directly 
in the hyperradial equations (|^). 

The strength and the range parameters 5*3 and 63 are now adjusted to repro- 
duce the position of the level at 0.38 MeV and in addition to get the ground 
state energy as close as possible to the observed value of —7.27 MeV. This re- 
quirement uniquely specifies the range and the strength of the three-body force. 
The resulting parameters are given in table 1. 

Properties of the system. The model is now completely defined and we can 
compute various measured properties of these two 0"^ states. The crucial quantity 
is the effective potential in the hyperradial equation, i.e. Wn = ^^/(2m)((A„ + 
^)/ p^ — Qnn) which is shown in fig. 1. The position of the level is indicated by 
the dashed line. The combined effect of the three-body centrifugal and Coulomb 
forces creates a strong barrier extending approximately from 13 to 60 fm. This 
barrier is responsible for the extremely small width of the resonance. 
The root mean square radius is in the three-body model given by 



Rrms = ^jnl + , {P') = E f fnip)p'dp , (7) 

where Ra = 1.47 fm is the root mean square radius of the a-particle. As seen in 
table 1 we underestimate the radius slightly in contrast to a larger overestimate 
found in other three-body calculations. This difference is due to the attractive 
three-body force which tends to decrease the size of the system. 
The monopole matrix element is written as 

3 O O POO 

M = 2E(0+ |r,2|0^) = ^(0+|p^|0+) = i-E / fi''\p)fi'^\p)p'dp , (8) 

i=l n •'^ 

where f^^^ (p) and f^^"* (p) are the radial wavefunctions for the ground state and 
the excited state, respectively. The numerical result given in table 1 is found with 
the wavefunction for the excited state 02" calculated as a true bound state by using 
an infinite wall at 30 fm. We overestimate the matrix element significantly less 
than other calculations indicating that we obtain a rather good description of the 
spatial configurations of the two 0"*" levels. 

To estimate the width of the 0^ state we first use the quasiclassical ap- 
proximation, that is calculate the penetrability of the combined centrifugal and 
Coulomb barrier in the lowest adiabatical channel. The penetration factor (~ 
10~^) is multiplied by the number of assaults per unit time on the barrier (~ 
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1 MeV/%) computed as the average radial velocity divided by two times the 
width of the potential well. In this approximation the width is about 100 eV. It 
is already extremely narrow compared to the resonance energy of 0.38 MeV, al- 
though still an order of magnitute larger than the experimental value. The main 
source of the underestimate is in the neglect of the higher adiabatical channels, 
especially the second one which has an attractive pocket at p ~10 fm, see fig.l. 

We now find the complex energy solution of the system (^) with the boundary 
condition (^. In order to obtain a converged result we have to integrate the 
radial equations numerically up to p ~80 fm. As we can solve numerically the 
Faddeev equations (|l|) only until p ~40 fm, we have to extrapolate the eigenvalues 
accoriding to equation @). Because of the low-lying ^Be resonance the uncertainty 
in the third term of the asymptotic expansion leads to an uncertainty in the width 
of about 10 eV. Within this uncertainty the resulting width is r=20 eV which is 
remarkably close to the experimental value. 

The width of the 0^ level in ^^C has also been considered using a two-body 
^Be+a approximation within the a-cluster model (see, for example, [jrB|). In 



this approximation the ^^C first decays into '^Be and an a-particle, then ^Be in 
turn decays into two a-particles. In our genuine three-body approach the decay 
proceeds directly into the three-body continuum and the narrow resonance ^Be 
comes into play implicitly through the structure of the angular eigenfunctions 
and eigenvalues. 

The spatial shapes of the states are somewhat difficult to visualize due to the 
symmetry of the three particles. The wavefunction depends on six coordinates, 
which conveniently now can be chosen as the distance between two of the particles, 
the distance between their center of mass and the third particle, the angle between 
these two directions and three remaining angles defining the orientation of the 
total system. The latter three are rotational degrees of freedom and they define 
a plane containing the three a particles. The wavefunction of a 0"*" state is 
independent of these rotationsl coordinates. It depends upon the first three 
variables, which define a triangle with a particles in the corners. To visualize 
the shape of the state we sample a sequence of these triangles from the density 
distribution. That is in this sequence the number of triangles of a certain shape is 
proportional to the absolute square of the wavefunction of the set of coordinates 
corresponding to this triangle. 

We then calculate the positions of the alpha particles for each set of coordi- 
nates in the sequence and mark them as points on a plot. The distribution of 
these points refiects the spatial structure of the corresponding state. The frame 
in which the positions of alpha particles are plotted is chosen in the following 
way: i) the center of mass of the system is positioned in the center of the plot; ii) 
the system is oriented so that two of the three particles are in the lower halfplane 
and the last particle is in the upper halfplane; iii) one of the two principal axes 
of inertia is aligned vertically. Note that the condition (ii) introduces certain 
up-down asymmetry. For symplicity the origin of the system of coordinates on 
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the plot is put at the lower left comer. 

The result is shown in fig. 2 for the ground state, which clearly resembles 
an equal side triangle - a consequence of the fact that the ground state contains 
predominantly the lowest and simplest angular eigenfunction. The corresponding 
contour plot for the excited state is shown in fig. 3. Unlike the ground state the 
excited state has larger contributions from the higher angular eigenf unctions. The 
structure is therefore a complicated interplay between different configurations 
with the tendency to add up to a somewhat elongated triangle. 

In addition we calculate the maximum probability configuration. The latter 
for the ground state is close to an equal side triangle with the length of the sides 
equal to 2.98 fm. For the excited state it is a somewhat elongated triangle with 
the length of the sides equal to 3.9, 4.2, 6.1 fm. 

Conclusion. We have formulated an approach to deal with the three-body 
Coulomb problem in the continuum. The approach is then applied to three 
interacting a-particles which gives a description of the lowest two O"*" levels in 
^^C. The second of these states is in the 3a-particle continuum held together by 
the Coulomb barrier. We construct adiabatic hyperradial potentials and estimate 
the width of the excited state by solving the system of hyperradial equations for 
the complex energy. The width, estimated for the first time consistently within 
the three-body model, is found to be 20 eV which is in good agreement with the 
experimental value proving the validity of the formulated approach. The model 
includes a genuine three-body force, which imitates off-shell effects or residual 
three-body interactions. Furthermore, the model constitutes an improvement 
over previous three-body calculations when we compare with measured values of 
the root mean square radius and the monopole transition matrix element. The 
spatial structure of the ground state is confirmed to be predominantly an equal 
side triangle while the excited state has a much broader probability distribution 
with a tendency to look like a fiat triangle. 
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Figure Captions 

Figure 1 The effective potentials Wn = {{Xn + 15/4)/p^ — Qnn)^^ /{2m) for the 
three lowest adiabatic channels as functions of hyperradius. The dashed 
line represents the position of the O2" resonance. 

Figure 2 Contour density plot of the ground state described in detail in the text. 
The unit on the axes is fm. 

Figure 3 The same as in fig. 2 for the excited state. 
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Table 1 

Parameters and properties of the 0"*" states. The resuhs of this work are given in 
the first row. The results from the calculations in ref. 



Ill and ref. 12 



are given 

in the next two rows and the measured values are shown in the last row. We 
show the parameters (strength 5*3 and range 63) of the three-body interactions, 
which lead to the energies E{Of) of the two O"*" states measured relative to the 
threshold of three free a-particles. Finally we show the root mean square radius 
Rrms of the ground state, the monopole matrix element M and the width F of 
the excited state. 
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